High-speed observation of the piston effect near the gas-liquid critical point 
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We measured high-speed sound propagation in a near-critical fluid using a ultra-sensitive inter- 
ferometer to investigate adiabatic changes of fluids on acoustic timescales. A sound emitted by very 
weak continuous heating caused a stepwise adiabatic change at its front with a density change of 
order 10 -7 g/cm 3 and a temperature change of order 10 _5 deg. Very small heat inputs at a heater 
produced short acoustic pulses with width of order 10/xsec, which were broadened as they moved 
through the cell and encountered with the boundaries. The pulse broadening became enhanced near 
the critical point. We also examined theoretically how sounds are emitted from a heater and how 
applied heat is transformed into mechanical work. Our predictions well agree with our data. 
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Thermal equilibration in one-component fluids takes 
place increasingly faster near the gas-liquid critical point 

at fixed volume [JESSSIa0|Il![l3, des P ite 
the fact that the thermal diffusion constant D tends to 
zero at the criticality. This is because the thermal diffu- 
sion layer at the boundary expands and sounds emitted 
cause adiabatic compression and heating in the whole 
cell after many traversals in the container. This heating 
mechanism is much intensified near the critical point due 
to the critical enhancement of thermal expansion of the 
layer. If the boundary temperature T w is fixed, the inte- 
rior temperature approaches T w on the timescale of the 
piston time 0, 



h = L 2 /4( 7 - 1) 2 D, 



(1) 



where L is the cell length and 7 = C p /Gv is the specific- 
heat ratio growing near the critical point. This time is 
much shorter than the isobaric equilibration time L 2 /AD 
by the very small factor (7 — l) -2 [Tl|. 

The previous experiments have detected only slow tem- 
perature and density changes in the interior region on 
timescales of order 1 sec. The aim of this letter is to 
report ultra-sensitive, high-speed observation of sound 
propagation through a cell filled with CO2 on the criti- 
cal isochore close to the critical point T c = 304. 12K. We 
can detect density changes of order 10 _8 g/cm 3 taking 
place on a timescale of 1/xsec. In confined fluids, adia- 
batic changes are usually caused by a mechanical piston, 
but they can also be achieved by a heat input through 
a boundary 0, H H, H S S 0. The second thermo- 
acoustic method is particularly efficient near the critical 
point. These adiabatic processes are very fundamental 
and ubiquitous, but we are not aware of any experiments 
detecting the underlying fast acoustic processes. 

Our experimental set-up is displayed in Fig.l, where 
the fluid temperature was controlled with precision of 
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FIG. 1: Front view of experimental set-up. Cell length L 
is either 10.3mm or 5.5mm. A dc current is sent to a film 
heater placed 0.5mm below the top boundary. Light pass- 
ing through the slit of width 0.25mm is used to detect small 
density changes using interferometry. Gravity is downward. 



±lmK using a Pt-thermometer and a four-wire resis- 
tance bridge |l2j . The upper and lower plates have area 
1 cm 2 and are made of Cu with high thermal conduc- 
tivity Xcu/ks = 2.8 x 10 23 /cm sec. The cell length L 
was either 10.3mm or 5.5mm. A thin NiCr-foil heater 
with thickness 3/im was placed d = 0.5mm below the 
top plate. The heat capacity of the foil is so small 
(Ch/kB = 8.4 x 10 19 /cm 2 per unit area) and the gen- 
erated heat was almost released to the fluid [13|, The 
side walls are made of stainless steel, whose thermal con- 
ductivity is 3.9% of that of Cu. We thus neglect heat 
flow to the side walls. A laser beam was sent at the cell 
center and a Twymann-Green interferometer was used 
to detect small density changes in a slit region of width 
0.25mm. Individual signals were very noisy and the data 
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points in the following figures are the averages over 418 
heat pulses successively generated in 2 seconds. Using 
sound pulses we first measured the sound speed c. For 
both L — 10.3mm and 5.5mm, our data on c closely agree 
with previous ones in the range T/T c — 1 >10 -4 [14(, but 
become independent of T — T c closer to the critical point 
because of the frequency-effect (inherent to short-time 
pulses) and/or the gravity effect 15]. 




heating 



FIG. 2: Normalized density change 8p(t)/p at the cell center 
in the time region < t < 0.4msec for T — T c — 150mK and 
30mK, produced by continuous heating in < t < 0.2msec in 
a cell of L = 1.03cm. Inset: Long-time behavior for T — T c — 
150mK in the time region 0.2 < t < 1.4msec. 
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FIG. 3: Profiles of pressure deviations propagating at the 
sound speed at T — T c = 150mK. Heating is stopped at t = 
200/isec, which results in a trapezoidal pulse (right). 

In Fig. 2 we show acoustic density variations Sp(t) at 
the cell center x — L/2 in a cell of L — 1.03cm at 
T-T c = 150 and 30mK, where c ^ 140 and 123m/sec, 
respectively. Continuous heating was applied in the time 
region < t < 0.2mscc and a sound signal arrived at 
the cell center only for t > (L/2 — d)/c ~ 40/^sec. The 
total heat supplied was Q = 367erg. The inset of Fig. 2 
displays a rapid relaxation at T — T c = 150mK after 
switching-off the heater. Fig. 3 illustrates spatial profiles 
of the pressure deviation 5p, whose inhomogeneity re- 
mains small even in the thermal diffusion layers. No- 
tice that dp takes large negative values in the layers. 
The thickness of the layers was only 2.6 x 10~ 5 cm at 
t = 0.2msec for T ~T C — 150mK. Right after switching- 



on the heater, step-like sounds were emitted on both sides 
and the one moving upward was reflected at the top at 
t = d/c. The two pulses merged for t > 2d/c ~ 10/xsec. 
After sweeping of the sound, we can see a stepwise in- 
crease in the density at t = L/2c,3L/2c, and 5L/2c. 
After the switching-off at t = 0.2msec, a trapezoidal 
pulse appeared moving back and forth in the cell. The 
trapezoidal part should be diminished if the heating time 
would be a multiple of 2L/c. The interior temperature 
was slightly increased by ST — (dT/dp) s Sp and a heat 
flow was induced through the boundaries. In the present 
case ST ~ 0.2TSp/p ~ 10 _5 deg. As can be seen in the 
inset of Fig. 2, the subsequent relaxation in the interior 
region was very rapid. It is the same as that after cool- 
ing of the boundaries on long times scales (t^> L /a) , so 
it is an adiabatic process caused by contraction of the 
thermal diffusion layers. 

In Fig. 4 we show Sp(t) after very short heating in a 
cell of L = 5.5mm at T - T c = 500 and lOOmK, where 
c = 149 and 133m/sec, respectively. A dc electric cur- 
rent of duration time 4.5/isec passed through the heater, 
where the pulse shape was determined by the relaxation 
time 2/isec of the current amplifier. The total heat sup- 
plied was Q = 129erg. The pulse directly leaving down- 
ward and that reflected at the top can be distinguished 
since their peaks are separated by 2d/c (which is 6.7/isec 
for T — T c = 500mK). Afterward, the two-peak pulse thus 
formed moved with speed c in the cell and was gradually 
flattened. In Fig. 5 we show pulse propagation and its de- 
cay at T — T c = 500 and 35mK in the cell of 1.03cm on a 
longer time scale, where the pulses are singly-peaked be- 
cause of the longer pulse-duration time 7/xsec (> 2d/c). 
The pulse broadening and the tendency of homogeniza- 
tioa become more enhanced on approaching the critical 
nt. In these experiments, the fluid velocity in the 
se region was of order v = cSp/p ~ 10 _2 cm/sec and 
fluid displacement was of order Ax = vAt ~ 10 _7 cm 

For 



10/isec being the pulse-duration time. 



po: 
-pri 

the 

wilh At 

"such extremely small displacements the film heater did 
not disturb the sound propagation [Hil ]. We are inter- 
ested in the adiabatically increased energy 2? a d(= pAx) 
in the pulse region per unit area. To linear order it is 
expressed as 



P f PC f 

-E a d = — dxSp(x, t) — — / dtSpi 

P J P J 



x,t). 



(2) 



The ratio of E^d to the total heat supplied Q — J dtQ(t) 
represents the efficiency of transforming applied heat to 
mechanical work. In our experiments E at \/Q was about 
0.11 — 0.12 at the first arrival of pulses at the cell center 
over wide ranges of T — T c and Q. However, a pulse has a 
tail persisting in later times and there remains ambiguity 
of order 10% in the calculation of i? a( j. 

We give theoretical interpretations of the above results 
neglecting the effect of the bulk viscosity. Let a small 
amount of heat be supplied to a near-critical fluid at a 
rate Q(t) per unit area for t > in a one-dimensional 
geometry. The volume expansion rate per unit area near 
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(4) 
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FIG. 4: 6p(t)/p at the cell center for T - T c = 500mK and 
lOOmK in the time region < t < 0.12msec, produced by 
pulse heating of width 4.5/isec in a cell of L = 5.5mm. 



■vfhere the Maxwell relation (dp^ 1 / ds) p — (dT/dp) s 
been used. The adiabatic coefficient (dT/dp) s is 
rly equal to the derivative along the coexistence curve 
(3T/<9p) cx near the critical point ^(j- Eor C0 2 it is 
to T c /6.98p c ^^ and the above relation becomes 
•out/p— 1.38xlCT 13 (5 with Q in erg/cm 2 sec. The heat- 
g rate used for our data in Fig. 2 was Q = 0.183 x 10 7 
then Sp ou t/p = 2.55 x 10~ 7 . This theoretical value 
irly agrees with the observed height of the first step 
2.2 x 10~ 7 in Fig. 2. The efficiency of energy transfor- 
ation discussed below Eq.(2) becomes 



a) 500 mK 

b) 35 mK 



FIG. 5: 6p(t)/p at the cell center for T - T c = 500mK and 
35mK in the time region < t < 0.6msec, produced by pulse 
heating of width 7psec in a cell of L = 1.03cm. 



the heater is 



Vi = 



dp- 1 } Q(t) 
ds 



T 



(3) 



where s is the entropy per unit mass. If the heater 
is attached to the upper plate, this volume change 
then produces a sound wave propagating downward as 
$Pcmt(t ~ x/c), where x is the distance from the heater 
and Sp out (t) is the density increase in the sound emitted 
at the heater at time t. Since acoustic disturbances move 
with the sound speed c, the mass conservation gives the 
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pfdT 
T\dp 



(5) 




which is 1/6.98 = 0.14 for near-critical C02- On the 
other hand, the experimental values were in the range 
0.11 - 0.12. 

As a complicating factor in our experiment, however, 
the heater is placed d = 0.5mm below the top plate. If 
V> is the volume expansion rate in the thermal diffusion 
kyer at the top plate, we require cSp on t(t)/p = Vi + V%. 
Ir. the time region where {Dt) 1 / 2 <C d <C ct < L, some 
calculations yield a modified convolution-type formula, 



Spout(t) 



cT 



dT 
dp 



dt'Q(t-t')iP(t'/t 2 )/t 2 , (6) 



w|here t 2 = L 2 /4c 2 ii is a characteristic time of heat ex- 
ange between the boundary and sounds in the range 
c <C t% <C t\. The scaling function ip(u) is positive 
and its integral ip(u) — J " dvip(v) tends to 1 for large u. 
TJhe function ip(u) itself appeared in the original theory 
? Thus Eq.(5) reduces to Eq.(3) when Q(t) changes 
much slower than t 2 - 

We clarify the relationship of the formula (3) and the 
assumption in the theory of the piston effect |2( . Let the 
heat input Q(t) changes slowly compared with the acous- 
tic time t a — L/c. We suppose a time interval with width 
St t a , in which Q(t) is almost unchanged. After many 
sound traversals, the adiabatic pressure and density in- 
creases in the interior region are given by 



2 2 6t p (dT\ . 



(7) 



as a result of superposition of many steps. In terms of 
the incremental heat supply SQ — QSt we find 



SQ 
nTV 



(8) 



where the Maxwell relation (dp/ds) p = p 2 (dT/dp) s has 
been used. If 5p is assumed to be homogeneous, Eq.(8) 
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follows from the space average of 5p = (dp/ds) p Ss + 
(dp/dp) s Sp 0], since the space integral of Ss is SQ/pT 
and that of 5p vanishes. In addition, the average tem- 
perature deviation is written as (5T) = SQ/CyL where 
Cy = pT{ds/dT) p is the isochroic specific heat [To| . 
From Eq.(3) we may thus reproduce the original theory 
of the piston effect. 

In summary, we have first measured rapid acoustic 
responses in near-critical CO2 to extremely small heat 
inputs, which constitute a basis of understanding adi- 
abatic processes in fluids. We demonstrate that adia- 
batic changes take place on the timescale of L/c. As 
Fig. 2 suggests, it is intriguing that a homogeneous adi- 
abatic change can be achieved almost instantaneously if 
the heating time is a multiple of 2L/c. The theoretical 
expressions (3) and (4) are useful to analyze these results 
and can reproduce the original theory of the piston ef- 
fect valid on long timescales. The theoretical height of a 
stepwise sound in Eq.(4) and efficiency E ad /Q in Eq.(5) 
fairly agree with the experimental values. Differences of 
20 — 30% remain, probably because a fraction of heat 
escaped to the stainless part and the pulse shape was 
considerably broadened even at the first arrival at the 
cell center. 

There are a number of unsolved problems, (i) We 
should understand how a sound pulse is reflected at a 
boundary wall and how it is damped in the bulk re- 
gion. In our case, the change of the pulse shape on 



reflection is apparently more marked than that due to 
the bulk damping, as can be inferred from comparison 
of the results for L = 5.5mm and those for L = 1.03cm. 
This is still the case even very close to the critical point, 
(ii) The bulk viscosity exhibits strong critical anomaly 
as C — 0.03pc 2 ^ 0| ) where t% is the order parameter 
relaxation time (18/xsec at T - T c = 30mK) [H). It 
gives a bulk acoustic damping, but it also affects the 
thickness of the thermal diffusion layer when the typ- 
ical timescale is longer than j(/pc 2 = 0.037t{ @, [§]• 
The effects of the bulk viscosity on heat transport are 
not yet well investigated, however, (iii) Upon heat ex- 
change the wall temperature is not fixed when the effu- 
sivity ratio a w = (X W C W j 'ACp) 1 / 2 is small 0, 0- Here 
Am, A, C w , and C p are the thermal conductivities and 
the specific heats (per unit volume) of the solid and 
fluid, respectively. For finite a w the piston time is 
given by h' = [(1 + 0/(7 - 1)} 2 L 2 /4D 3], which 
becomes t\ = C 2 -L 2 /AX W C W for a w <C 1. In our case 
a w = 3 x 10 3 (T/T C - l) 92 between Cu and C0 2 and 
a w < 1 is reached for T — T c < 50mK. This crossover 
should be studied. 
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